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Abstract 

We carry out a detailed investigation of neutral {D^) and charged {D~) impurity states of 
hydrogen-like donors in spherical semiconductor quantum dots. The investigation is carried out 
within the effective mass theory (EMT). We take recourse to local density approximation (LDA) 
and the Harbola-Sahni (HS) schemes for treating many-body effects. We experiment with a variety 
of confining potentials: square, harmonic and triangular. We observe that the donor level undergoes 
shallow to deep transition as the dot radius {R) is reduced. On further reduction of the dot radius 
it becomes shallow again. We term this non-monotonic behaviour SHADES. This suggests the 
possibility of carrier "freeze out" for both and D~ . Further, our study of the optical gaps 
also reveals a SHADES transition. 
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I. INTRODUCTION 



Impurities play a central role in semiconductor technology. Performance of a semicon- 
ductor device is dictated by shallow dopants. Shallow dopants like phosphorous and boron 
alter the conductivity of bulk silicon by several orders of magnitude. Deep defects on the 
other hand are known to degrade device performance. An interesting point is that the con- 
centration of these impurities is miniscule: less than a tenth of a percent. Relevance of 
the role of shallow donors in semiconductor quantum dots (QDs), which are essentially zero 
dimensional structures, thus cannot be overestimated. 

Nearly two decades ago Bastard [ij reported the first calculation of the binding energy 
of a hydrogenic impurity in two-dimensional quantum wells (QW). Bryant extensively 
studied hydrogenic impurity states in one-dimensional quantum well wires (QWW). Zhu et 
al. j^l were the first to study the effect of hydrogenic impurity in spherical quantum dots 
(QDs). The system chosen was {GaAs — Gai_xAlxAs). Calculations based on the hydrogenic 
impuritvstates in presence of electric or magnetic field have also been reported over the past 
decade p, dill- Ranjan and Singh ^ reported studies of neutral donor states in II- 
VI and III-V QDs. They examined the effects of a size dependent dielectric constant e{R), 
where R is the radius of the spherical QD. As long back as 1990 Pang and Louie P| discussed 
the negative donor centers in semiconductor quantum wells. They calculated the binding 
energy of negative donor state {D~) in presence of magnetic field, solving the effective mass 
model exactly bv diffusion quantum Monte Carlo method. 

Zhu et al. |lol . \v\ worked out the effect of dimensionality and potential shape on the 
binding energy of and states of the impurity. Further, they also studied the effect of 
electron correlation in the binding energy. They showed that the shape of the confinement 



is less important than the confinement strength. They also mentioned the importance of 
electron correlation in their results. Ranjan and Singh ^ made a systematic study of the 
non-monotonic shallow deep transition (SHADES) of the impurity in a semiconductor QD. 
SHADES occurs when one decreases the size of the QD. First the impurity level becomes 
deep and then below a critical size it becomes shallow again. 

In the present work we study both the neutral D° and the negatively charged D~ hydro- 
genic donor states in a three-dimensional spherical quantum dot (QD). We study the effect 
of the shape of confinement potential on the binding energy of these donors. Many-body 
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effects are taken into account by employing the density functional methodology. To treat 
the exchange effects accurately, particularly for negative ions, we employ the Harbola-Sahni 

compare the HS results 

with those of the local density approximation (LDA). 



II. MODEL 



Our calculations have been performed on a spherical QD We have solved the Schroedinger 
equation self-consistently within the effective mass theory (EMT) approximation. We have 
included electron-electron interaction effect within the local density functional approach 
and the Harbola-Sahni scheme. The Hamiltonian of the system is given by the following 
expression 



]\f N N 
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The first term in the above Hamiltonian is the kinetic energy of the electrons, the second 
term is the electron-electron interaction energy, the third term is the external potential for 
N-electron system and the last term is the impurity potential. This Hamiltonian is written 
in Hartree's atomic unit, where h = rrie = e = 1 and unit of energy is 27.2 eV and that of 
length (size) is 0.53 A. In the above equation m* is the effective mass of the electron inside 
the QD in units of rrie, the free electron mass. The above many-body Hamiltonian can be 
reduced to a single particle Schroedinger type Hamiltonian for the ground state within the 
Kohn-Sham formulation of density functional theory (DFT) |l3|. 

In order to study the shape dependence of binding energy, we model 
potential as follows: 

/ (Vo/R'y -Vo r<R 
Vext{r) = < (2) 
1^ r> R 

where Vq is the depth of the potential. This can be given by the conduction band offset 
(valence band offset) between the QD and the surrounding layer for the electron (hole). R is 
the radius of the QD and k assumes positive integral values from 1 to a very large number. 
Changing the value of k results in the change of the shape of the potential. In particular 
k = 1 is quasi-triangular, A; = 2 is quasi-harmonic confinement and > 10 is quasi-square 
well confinement. As /c oo, the potential becomes a square well. 
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FIG. 1: A schematic diagram, shows the impurity doped spherical quantum dot (QD). The left 
side of the figure depicts a spherical QD of size R, surrounded by a dielectric (glasses, polymers, 
organic solvents or oxides and hydrides). There is a hydrogenic impurity at the center of the QD. In 
our model we vary the shape of the confinement, namely quasi-triangular {k = 1), quasi-harmonic 
{k = 2) and quasi-square {k > 10) (see Eq. ((21) )• The right hand side of the figure shows one typical 
model potential, corresponding to the shape index /c = 80 (quasi-square well). The depth of the 
well is Vq and the diameter is 2R. The neutral and negative donor levels are schematically shown 
along with the conduction band minimum (CBM) with corresponding energies E{D^), E{D~) and 
£"0 respectively. Their respective binding energies Eb{D^) and Eb{D~) are indicated by arrows. 

The total ground state energy of an N-electron system in terms of the density is given as 

E [p] = Ts [p] + Eh [p] + E,, [p] + E,,t [p] + E^^^ [p] (3) 

where Ts[p] is the kinetic energy functional of a non- interacting many electron system, 
Eh[p\ is the Hartree energy, E^cIp] is the exchange-correlation energy functional, E^xtyp] 
is the energy functional due to the external potential and Eimp[p] is the energy functional 
because of the impurity potential. Note that p(r) is the charge density. As pointed out in 
Sec. Ewe have used two different approaches for the calculation of the exchange-correlation 
potential (\4c), namely the LDA in its Gunnarsson-Lundqvist ^| parameterized form and 
the exchange potential using HS scheme. Within the HS scheme, the exchange potential is 
interpreted as the work done in moving an electron in the field of its Fermi hole. 



The HS scheme calculates nearly exact exchange interaction. For the atomic case, the 
method essentially reproduces Hartree-Fock results for all the elements in the periodic table 
[isl l. Also, HS scheme has self-interaction correction built in it and hence calculates very 
good eigenvalues and total energies for a single electron system, whereas the LDA is known to 
yield poor results for the hydrogen-like atoms due to non-cancellation of the self-interaction 
of the electron. 

A hydrogenic impurity in a host bulk semiconductor gives rise to a shallow donor state 
in the energy gap of the host. We refer to this as neutral donor state (D^) and its energy 
as E{D^). The energy required to take the electron from the neutral donor state into the 
conduction band is termed as binding energy, Eb{D^) and is of the order of 10 - 50 meV. 
Hence the binding energy of a neutral donor can be written as 

Eb{D'')=Eo~E{D') (4) 

where Eq is the ground state energy of an electron in the conduction band. 

A negative donor center in a semiconductor is formed by a neutral donor center trapping 
an extra electron. The binding energy of a negative donor can be interpreted as follows: We 
take an electron from D~ to the conduction band. Initially the energy of the system was 
E{D^), the final energy is Eq + E{D^). Therefore the binding energy of a negative donor 
can be given by subtracting the final energy from the initial energy, i.e. 

Eb{D-) = E{D') +Eo- E{D-) (5) 

where E[D^) is the ground state energy of a neutral donor and E{D~) is the ground state 
energy of a negative donor. This is the energy required to promote an electron from to 
the conduction band. Naturally the binding energy of a negative donor is smaller than the 
neutral donor. 

As we have mentioned earlier the HS scheme has self-interaction correction built-in, hence 
it calculates very good eigenvalues for small number of electrons. In fact the total energy of 
a single electron with and without the impurity potential is exactly equal to the eigenvalue 
of an electron with and without the impurity potential respectively. In other words, Eq = eo 
and E{D^) = ei, where eo is the eigenvalue of an electron in the conduction band minimum 
and ei is the eigenvalue of the neutral donor. The expression for binding energy of a neutral 
donor in terms of the eigenvalues is given as 

= eo-ei (6) 



Thus the binding energy of an electron in this case is also equal to the optical transition 
energy. We can also write the binding energy of a negative donor in terms of the eigenvalues. 
The total energy of a hydrogen-like negative donor is given by the formula 

1 f p{r)p{r') 



E{D- 



2t2 

2 



-drdr' + E^c[p\ - I V^c{r)p{r)dr 



(7) 



where €2 is the eigenvalue of two-electron system. From the definition of the binding energy 
of a negative donor (Eq. 

p{f)p{r') 



ei + eo - 2e2 



-dfdr' - E^c[p] + / V^c{r)p{r)df 



ei + eo - 2e2 + f/ - 2E^,[p] + j V^,{r)p{r)d7 



(9) 



where 



U 



p{r)p{r') 



dfdr' + E^c[p\ 



(10) 

is the effective electron-electron interaction energy. We call it the "Hubbard t/" [M Q- 
The last two terms in the Eq. © cancel each other (true only for one orbital system) in 
HS scheme but they do not cancel within LDA. Hence, within the HS scheme, the binding 
energy of a negative donor in terms of eigenvalues can be written as 



Eb{D- 



ei + Co — 2e2 + U 



(11) 



Within the LDA, U = Ucoui + Ux + Ucorr whereas within the exchange only HS scheme, 
U = Ucoui + Ux, where Ucoui is the Coulomb energy, Ux is the exchange energy and U corr is 
the correlation energy. 

III. RESULTS AND DISCUSSION 



A. The Binding Energy 

We present the results of binding energy of a hydrogen-like impurity [Z = 1). As men- 
tioned earlier we are employing Hartree units, in which h = rrie = e = 1 and unit of energy is 



6 




FIG. 2: The binding energy of a neutral donor {Eb{D^)) is plotted as a function of the size R of 
the QD. The shape chosen is quasi-square well {k = 80). The depths of the potential well are taken 
to be Vo = 80, 40, 20, 10, 5 and 1. The figure shows a non-monotonic shallow to deep (SHADES) 
transition of the binding energy as the dot size is reduced. In the inset we plot the conduction 
band edge or LUMO (solid line) and the impurity level (dashed line) as a function of size for Vq = 
10 and k = 80. Hartree units are employed. 

27.2 eV and that of length (size) is 0.53 A. Unless otherwise stated all the results presented 
are based on HS scheme. But the HS scheme within the one electron system {D^) is no dif- 
ferent from an electron moving in an impurity potential; essentially the Ranjan-Singh (RS) 
^. However RS scheme is perturbative and therefore its results are approximate. On the 
other hand, our results are numerically exact. We solve the one - dimensional Schroedinger 
equation in spherical co-ordinates. 

In Fig. 121 we depict the binding energy {Eb{D^)) of a neutral donor {D^) as a function 
of the size R of the QD for various different depths of confinement potential, namely Vq = 
80, 40, 20, 10, 5 and 1. The shape chosen is quasi-square well {k = 80). The plots indicate 
a monotonic increase in the binding energy as the size is reduced. As we decrease the size 
further, the binding energy reaches a maximum and then it decreases after a critical size. 
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Thus a non-monotonic shallow- deep (SHADES) transition of the donor level is observed. A 
qualitative explanation of SHADES transition is as follows: The binding energy is the dif- 
ference between the energy of conduction band minimum (CBM) and the impurity energy. 
As the dot size is decreased the CBM rises monotonically as 1/R'', where 7 G [1.2 : 2.0] 
[isj ]. On the other hand, initially the impurity energy remains relatively constant in value 
as R is decreased and then it increases rapidly. This is because the impurity wavefunction 
is localized and consequently the impurity charge density is highly confined. This charge 
density does not sense the boundary initially. We have confirmed this picture by a detailed 
examination of the energies and the wavefunctions. A typical plot of CBM (Eq) and the 
impurity ground state energy {E{D^)) is shown in the inset of Fig.|21 Even a visual examina- 
tion reveals that the difference between the two energies is non-monotonic with size. Below 
a critical size, which we define to be Rshades, the donor level becomes shallow again. This 
is because below Rshades, the kinetic energy expectation value increases much more than 
the potential energy expectation value. Note that for large R the binding energy approaches 
the bulk limit. 

The binding energy is large for large value of Vq. This is because of large confinement 
of the charge carrier. We also see that the critical size Rshades at which SHADES oc- 
curs increases monotonically with decrease in the depth of the well. These observations are 
recorded in Table HI We explain these as follows: For smaller depth of the confinement po- 
tential (for example Vq = 1), the strength of the confinement is less. Conventional textbooks 
define the strength of the potential as VqR^ . Because of this, the electronic wavefunction 
is more spread out. As one decreases the size, the wavefunction senses the boundary much 
before {i.e at larger size) than the one for which the strength of the confinement is large 
(for example Vq = 80). Hence for small Vq, the energy of donor level increases rapidly. This 
implies that the binding energy is small for small Vq and the SHADES transition occurs at 
larger size. An examination of Table IJ reveals that VoR^g^jt^j^^g is approximately constant. 
In other words, Rshades ~ Thus, the binding energy and its maximum depends on 

the well depth Vq. The latter is a surface related property depending on surface termination 
of dangling bonds, the dielectric coating, etc. Thus, with experience, it maybe possible 
to engineer the magnitude of the binding energy, making it shallow, intermediate or deep 
depending on one's convenience. It suggests the possibility of ''defect engineering" in QDs. 

We pause to note that the hydrogenic impurity in bulk semiconductor is quite shallow. 
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FIG. 3: The binding energy of a neutral donor (D^) is depicted as a function of the size of the QD. 
The depth of the well is chosen to be Vq = 10 and the shape of the confining potential is varied, 
corresponding to the shape index k = 80, 10, 4, 2 and 1. Hartree units are employed. 

The carriers in the shallow levels of hydrogenic impurities are easily promoted into the 
conduction band at room temperature, thus enhancing the conducting properties of the 
semiconductor by several orders of magnitude. The increase in the binding energy of the 
hydrogenic impurity with decreasing size indicates that carrier "freeze out" will occur. Thus 
a nominally shallow donor will "go deep" . 

In Fig. ini we explore the dependence of the binding energy on size as the shape of the 
confinement potential is varied from quasi-triangular {k = 1) to quasi-harmonic {k = 2) and 
to quasi-square {k = 80). The potential depth is kept fixed at Vq = 10. We observe that, 
in this case also, the binding energy increases monotonically as the size is reduced, reaches 
a maximum and then decreases. However, the critical dot size at which the maxima occur 
is relatively insensitive to the shape parameter (k). Further, the dependence of the binding 
energy maximum on the shape parameter {k) is non-monotonic. This is also displayed in 
Table ITTl We point out the possible significance of this result. When the charge distribution 
is uniform throughout the spherical QD, then the shape of the potential maybe taken to 
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be quasi-harmonic. If there is excess charge on the surface then the shape of the potential 
maybe taken to be quasi-square. Thus, depending on how the charges are distributed in a 
QD, it is possible to engineer the magnitude of the binding energy. Unlike in Fig. |2l where 
both the magnitude of binding energy maxima and Rshades are sensitive to the depth of 
the potential chosen, we see in Fig. |21 that they are relatively less sensitive to the shape of 
the potential. 

Notice from Fig. El and the Table |n] that the binding energy maxima occur at different 
sizes for = 80 {Rshades = 0.55) and A; = 10 {Rshades = 0.65). Here the magnitude 
of the binding energy maximum is larger for k = 10 {Eb{D^) = 2.63) than for k = 80 
{Eb{D^) = 2.56). Similarly, the binding energy maxima occur at slightly different sizes 
for A; = 10 {Rshades = 0.65) and k = 4 {Rshades = 0.70) but the magnitude of the 
binding energy maxima are roughly the same {Eb{D^) = 2.63 and 2.64 for A; = 10 and 
4 respectively). Also, notice that the binding energy maxima occur at larger sizes as we 
change the shape index k from 4 to 2 and then 1, however, the magnitude of the binding 
energy maxima decreases. 

Figure El also reveals that the binding energy is largest when quasi-triangular confinement 
is used. The binding energy of quasi-harmonic confinement is larger than the quasi-square 
confinement. This is because the charge carriers are more confined in quasi-triangular po- 
tential than quasi-harmonic and quasi-square potential. Similarly charge carriers are more 
confined in quasi-harmonic than the quasi-square potential. Thus one can tune the magni- 
tude of the binding energy by tuning the shape of the potential. 

Figure m shows the binding energy of a negative donor {D^) as a function of the size 
of the QD. The parmaeters used are the same as for Fig. El Once again a non-monotonic 
SHADES transition is observed. For very small size the binding energy becomes negative. 
This implies that the impurity level is resonant with the conduction band. Note, SHADES 
occurs at larger size when the depth of the potential is taken to be smaller. Finally at large 
size the binding energy approaches the bulk limit. 

If we compare the binding energy plot of (Fig. E|) with the binding energy plot of D~ 
(Fig.E}, we note the following: The binding energy of is smaller than for a fixed depth 
of the confinement Vq and shape index k. This is obvious because the D~ level is higher 
than the level by an amount approximately equal to U, the effective electron - electron 
interaction energy. For a fixed k and Vq, the SHADES transition for D~ occurs at larger size 
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TABLE I: The table shows the maximum binding energy of D and D~ states and the critical 
size Rs HADES at which the SHADES transition takes place with varied depth of the confinement 
potential (Vq). We note that over a large range of confining potential (Vq = 1 - 80), the strength of 
the confinement (Vq-R^) at Rshades remains relatively unchanged (column III and VI). The shape 
of the confinement potential chosen is quasi-square well {k = 80). Hartree units are employed. 
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FIG. 4: The binding energy of a negative donor {D~) is plotted as a function of the size of the QD. 
The shape chosen is quasi-square well {k = 80). The depths of of the potential well arc chosen to 
be Vb = 80, 40, 20, 10, 5 and 1. The figure shows a a non-monotonic SHADES transition. Hartree 
units are employed. 
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FIG. 5: The binding energy of a negative donor (D^) is depicted as a function of the size of the 
QD. The depth of the well is chosen to be Vq = 10 and the shape of the confining potential is 
varied, corresponding to the shape index k = 80, 10, 4, 2 and 1. Hartree units are used. 

than for D° i.e. Rshades{D^) > Rshades{D^)- This is because the wavefunction of the 
negative donor is more spread out due to the Coulomb repulsion whereas, the wavefunction 
of neutral donor is more confined. As the size becomes small, the wavefunction of negative 
donor senses the boundary relatively earlier and gets squeezed. As we decrease the size 
further, below a critical size {Rshades{D^)), the charge carriers are delocalized and hence 
leak out of the well. Therefore the binding energy of negative donor (D^) decreases at a 
size {RsHADEs{D~)) larger than the size {Rshades{D^)) at which the binding energy of 
neutral donor (1)°) decreases. 

We depict the dependence of the negative donor on the shape of the confinement potential 
in Fig. El The parameters used are the same as for Fig. IHl The binding energy exhibits a 
non-monotonic SHADES transition. The critical dot size at which the maxima occur is 
even more insensitive to the shape parameter (k) than it is for the neutral donor, as is 
evident from Fig. El and a comparison of columns II and IV of Table |nl The dependence 
of the binding energy maximum on the shape parameter (k) is non-monotonic as it is for 
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TABLE II: The table shows the maximum binding energy of D and and Rs hades at with 
varied shapes (k) of the confinement potential. The depth of the confinement potential chosen is 
Vo = 10. Hartree units are employed. 



_k RsHADEsjD'') Eb{D'') RsHADEsjD-) Eb{D^) 

80 0.55 2.562 0.80 0.536 

10 0.65 2.633 0.80 0.591 

4 0.70 2.645 0.90 0.638 

2 0.80 2.558 0.95 0.654 

1 1.05 2.329 1.15 0.619 



the neutral donor. The binding energy maxima occur at larger size as we change the shape 
index k from 2 to 1, however the magnitude of the binding energy decreases. These and 
other aspects are displayed in Table |n] Note the binding energy of the negative donor is 
much smaller than the binding energy of the neutral donor. 

We have also calculated the binding energies of neutral and negative donors using Eq. ^ 
and Eq. (jllj) . As pointed out earlier, we use two different methodologies, LDA and HS 
scheme. We have compared the binding energies of a neutral donor, obtained using Eqs. Q 
and We have also compared the binding energies of a negative donor, obtained using 
Eqs. (0) and (fTT|) . within the HS scheme. We find that the HS scheme gives identical values 
of Eb{D^) using Eqs. (jH) and (0). It also gives identical values of Eb{D^) using Eqs. (0) 
and (jllj) . We did the same exercise within the LDA. But LDA gives poor results for both 
cases. 



B. The Optical Gap 

Figure El shows the eigenvalue difference (eo — ci) between an electron in the conduction 
band minimum and the neutral donor. This difference (Ae) is representative of the optical 
gap 1^ of the system. The depth of the potential chosen is Vq = 10. Three different shapes 
of the confinement potential are chosen, corresponding to = 1, 2 and 80. We show results 
obtained using both the LDA and the HS scheme. We observe that the optical gap also shows 
non-monotonic SHADES transition. As noticed earlier, the optical gap obtained using the 
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FIG. 6: The eigenvalue difference of an electron in the conduction band minimum and the neutral 
donor (eo — ei) is plotted as a function of the size of QD. The depth of the potential is chosen to 
be Vo = 10. Three different shapes of the potential are chosen corresponding to A: = 1, 2 and 80. 
For comparison we plot results of the calculations of both LDA and the HS scheme. Hartree units 
are employed. 

HS scheme is identically equal to the binding energy of a neutral donor, Eb{D^) (Eq. (jlj)). 
Thus within the HS scheme, the optical gap follows the same track with the size (R) and 
the shape parameter k as shown for Eb{D^) in Fig. El The results obtained using LDA are 
different from those obtained using the HS scheme. The figure reveals that the optical gap 
calculated using LDA is smaller than that of the HS scheme. This is because magnitude of 
eigenvalue within the LDA for hydrogen-like systems is smaller due to the self-interaction 
of electron. 

Figure [3 shows the eigenvalue difference (eo — £2) between an electron in the conduction 
band minimum and the negative donor (optical gap). The parameters chosen are the same 
as for Fig. IHl For comparison we show results obtained using both the LDA and the HS 
scheme. Once again we observe the non-monotonic SHADES transition in the optical gap. 
The optical gap shows similar features with size and the shape of the confinement potential 
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FIG. 7: The eigenvalue difference of an electron in the conduction band minimum and the negative 
donor (eo — £2) is plotted as a function of the size of QD. The depth of the potential is chosen to 
be Vo = 10. Three different shapes of the potential are chosen corresponding to A: = 1, 2 and 80. 
For comparison we plot results of the calculations of both LDA and the HS scheme. Hartree units 
are employed. 

as those shown in Fig. |S| In contrast to the single electron case, we find that the optical gap 
calculated using LDA is larger than that of the HS scheme. 

In the context of the optical gap, we would like to elucidate the connection between 
our calculation and the work of Bhargava et a/.|l2^. The deep defect levels in QD would 
act as efficient traps for photo-excited carriers. Enhanced luminescence results when these 
carriers undergo radiative recombination. This scenario was proposed by Bhargava et al. 
who found that incorporating Mn impurities in ZnS nanocrystallites results in a spectacular 
enhancement of the luminescence efficiency (18%) with a corresponding decrease of lifetime 
from milliseconds to nanoseconds. Thus high luminescence is perhaps at least partly related 
to the deepening of impurity levels. The deepening of the nominally shallow hydrogenic level 
observed by us implies that we maybe able to use them as efficient luminescence center. 
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FIG. 8: This plot shows the dependence of the binding energy of negative donor (meV) on the QD 
size R {nm) and the potential well depth Vq. The shape of the potential chosen is quasi-square 
well {k = 80). The materials chosen are GaAS and ZnS. Note that for ZnS, the defect levels are 
very deep (0.2 - 0.3 eV). Hence carrier "freeze out" is a distinct possibility in ZnS. Also, note 
that the SHADES occurs at very small size {R < 1 nm ) in ZnS nanocrystals whereas, in GaAs 
nanocrystals SHADES occurs at i? ~ 10 nm. 

IV. CONCLUSION 



Shallow dopants play a critical role in semiconductor technology. The present work as 
encapsulated in Figs. El - El presents an interesting scenario. This scenario is presented for 
the technologically relevant case GaAs and ZnS QDs in Fig. |H1 Figure |H1 is for a negati vely 
charged donor where many-body effects are relevant. The inset is for neutral donor 21 1. 
However as noted earlier both oi D~ show similar behaviour. A nominally shallow donor 
{D^ or D^) becomes deep as the dot size is decreased. This implies that carriers will "freeze 
out". On further decrease of the dot size the dopants may once again become shallow 
(SHADES). An examination of Fig.jHlalso indicates that the binding energy and its maxima 
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TABLE III: This table shows the potential depths in eV for materials, namely GaAS and ZnS. 



Vo Vo (GaAs) in eV Vq (ZnS) in eV 

80 0.920 17.86 

40 0.460 8.93 

20 0.230 4.46 

10 0.115 2.23 

5 0.058 1.12 

1 0.012 0.22 



depends on the well depth Vq. The latter is a surface related property depending on surface 
termination of dangling bonds, the dielectric coating, etc. Thus, with experience, it maybe 
possible to engineer the magnitude of the binding energy, making it shallow, intermediate or 
deep depending on one's convenience. This suggests the so called "synthetic tailor ability" 



22j | of binding energy by selecting a suitable dot size or dielectric coating. This may usher 



in the possibility of "defect engineering" in QDs. As is clear from Fig. |H1 the doped GaAs 
QDs of size i? < 20 nm (200 A) will be susceptible to carrier "freeze out". Similarly the 
doped ZnS QDs of size R < 2 nm (20 A) will be susceptible to carrier "freeze out". 

We also observed the SHADES transition in the optical gap. We note the relevance of our 
calculations in Bhargava et a/.'s work [2^. In this work, they proposed that the deep levels in 
materials such as Mn doped ZnS would act as efficient traps for photo-excited carriers. This 
leads to enhanced luminescence when these carriers undergo radiative recombination. If we 
accept this point of view, then a suitably tailored "deep" dopant may act as a luminescent 
center. 

We note that we are working in EMT-LDA approximation. Preliminary work on an 
elaborate tight binding calculation on neutral donor seems to indicate SHADES behaviour 
2^. The donor impurities in the present work are located at the center of a spherical 
quantum dot. It would be interesting to carry out investigations in which the shape of the 
QD is varied and the donor is moved off center. It would also be of interest to investigate 
helium-like donors, e.g. S in Si. In this case many-body effects would be significant, since 
one would have at least three charge states of the donor, namely D~^, and D~ . We plan 
to undertake some of these tasks in future. 
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